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ABSTRACT 

We present proper motion of the sub-parsec scale jet in a nearby elliptical galaxy 3C 66B. 
Observations were made using the VLBA and partly Effelsburg 100-m telescope at 2.3 GHz and 
8.4 GHz at 10 epochs over 4 years. The 8.4 GHz images showed that a proper motion increases 
from 0.21 to 0.70 mas year~^, corresponding to an apparent speed of 0.30 c to 0.96 c, with 
a distance from the core on a sub-parsec scale. Our investigation suggests that the apparent 
increase of the proper motion can be explained by changes in the viewing angle, according to a 
relativistic beaming model. However, we still cannot eliminate the possibility that acceleration 
of the jet outflow speed or of changes of cmissivity profile in the two-zone jet might be found in 
3C 66B. 

Subject headings: galaxies : individual (3C 66B) - galaxies : active - galaxies : jets - galaxies : nuclei 



Introduction 

3C 66B is a giant elliptical galaxy with a red- 
shift (z) of 0.0213 (iHuchra. Voeelev. fc Gelle7 



vation at 5. GHz reveals the prese nce of a faint 



counterjet (jGiovannini et al. 1120011 ). The inten- 



1999t ) within the cluster A34 7 (iFanti et al. Ill982l) 



and t he Perseus super cluster (jGregorv. Thompson. 
19811 ). It has a prominent face-on dust lane, which 
i ndicates the occurrence of galaxy merger events 
(|Verdoes Kleiin et al. I [19991) . This source is also 



sity ratio between the jet and counterjet was es- 
timated to be approximately 10 within a distance 
Cjf 2 mas from the core, indicating possible accel- 



^^-liffyration of the jet outflow on a parsec scale. The 



known as an FR I radio galaxy and it has a strong 
jet and counterjet that extend to about 100 kilo- 
parsecs (kpc), which ha. ve been ob served at radio 
( Hardcastle et al. 1996), in frared (^ Tanslev et al. 
.20001 ), optical (iMacchetto et al. . 199ll) . and X-rav 
waves (|Hardcastle et al. 1120011) . 



orbital motion of the compact radio core in 3C 
66B, which has a period of 1.05 ± 0.03 years, 
was observed with a position accuracy of 10 micro 
ar cseconds (/Ltas) u sing phase-referencing VLBI 



( Sudou et al. 112003 ) 



High-resolution imaging using very long base- 
line interferometry (VLBI) has revealed an asym- 
metric core-jet struct ure similar to that ob- 
served in radio galaxies (iFraix-Burnet et al. Ill997 
Xu et al. I I2OOOI: iGiovannini et al. Il200ll ). Obser- 



Here, we report the detection of the proper mo- 
tion of the sub-parsec scale jet in 3C 66B based on 
multi-epoch observation over 4 years. Observa- 
tions and data reduction are described in §2, and 
measurements of proper motion, including the cal- 
ibration of the systematic shift of the position of 
the core are presented in §3. §4 focuses on the 
increase of the apparent velocity of the jet with 
distance from the core. 
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In this paper, we use a Hubble constant {Hq) of 
71 km s~i Mpc"i , a matter density (I^m) of 0.27, 
and a vacuum energy (r^A) of 0.73, which results 
in an angular size or separation of 1 milliarcsecond 
(mas) corresponding to 0.436 pc at the distance of 
3C 66B. 

2. Observations and Data Reduction 

3C 66B was observed at 2.3 GHz and 8.4 
GHz using VLBI during the period of 2001.20 
to 2002.45 (6 epochs) and of 2004.80 to 2005.69 (5 
epochs). In this paper, we name these two periods 
Period 1 and Period 2. Participating antennas for 
each epoch are listed in Table 1. We used the 
NRACQ VLB A (|Napier et al. | ll994') for the first 
10 epochs and the VLBA and Effclsbcrg 100-m 
telescope of the MPIfR for the last epoch. Due 
to bad weather conditions and antenna problems, 
the 8.4 GHz data for the second epoch (2001.48) 
showed very low quality in reproducing the faint 
jet structure, and so we removed this data from 
further imaging analysis. In addition, we only 
observed at 8.4 GHz for the last epoch (2005.69). 
Consequently, we obtained 10 sets of data for each 
frequency. 

The data were recorded in the VLBA for- 
mat with 2 intermediate frequency bands with 
a bandwidth of 8 MHz each for each frequency. 
The total bit rate for both frequencies was 128 
Mbps using a 2-bit sampling mode. The ob- 
served data were correlated with the VLBA cor- 
relator in Socorro. The correlated data for 3C 
66B were calibrated using the AIPS package de- 
veloped by NRAO. The data for all epochs were 
treated in a similar fashion; the residual phase 
delays and delay rates that were mainly due to 
antenna position errors and atmospheric phase 
fluctuat ions were calibrated usi ng global fringe- 
fltting (jSchwab &: Cotton I Il983f ) . Visibility am- 
plitudes were calibrated using measured system 
temperature and antenna gain curves. Bandpass 
calibration was performed using 0133-1-476 for all 
epochs. 

The calibrated visibility data were coherently 
averaged over 10 - 30 seconds according to data 
quality, which varied between epochs due to the 



different weather conditi ons. The imag es were 
made with DIFMAP (jShepherd I 119971 ) using 



^The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under coopera- 
tive agreement by Associated Universities, Inc. 



CLEAN and self-calibration. 

Figures [T] and [2] show the final images for 2.3 
GHz observation in Period 1 and 2, respectively, 
and Figures |3] and |4] show those for 8.4 GHz obser- 
vation in Period 1 and 2, respectively. Since the 
total observation time in the scheduled scan was 
almost the same for all the observations (about 
9 hours), the obtained dirty beams are similar to 
each other. However, dirty beams for the 2005.35 
epoch were not the same because of data loss at 
the St. Croix station due to antenna problems. 

3. Results 

3.1. Proper motion of the jet 

For quantitative evaluation of the source struc- 
ture and proper motion, we applied model fitting 
to several components for the visibility data using 
the MODELFIT program in DIFMAP. The core 
could be reasonably fitted with a single Gaussian 
component. On the other hand, since the jet did 
not show a discrete knot structure clearly in our 
images, Gaussian fitting of the jet region led to a 
single, largely extended component over the whole 
jet. In order to obtain the distribution of the in- 
dividual knot positions in the jet, we applied a 
point source fitting. The procedure for model fit- 
ting of the core and jet components is shown as 
below. First we carried out the Gaussian fitting of 
the core. In the residual map subtracted the fitted 
components, we found the local brightest peak in 
the jet region and measured its position and fiux 
as an initial model. Then we carried out the Gaus- 
sian fitting of the core and the point-source fitting 
of the jet components simultaneously. This pro- 
cess was repeated until the flux of the local peak 
became less than 3 a level. The xt values of the 
fitting results are calculated by using errors de- 
rived from the scatter of the visibility amplitude 
data within solution interval which is typically 30 
seconds. 

We successfully fitted 4 or 5 components in the 

2.3 GHz images and 2, 3, or 4 components in the 

8.4 GHz images. 

The errors for model fitting were estimated 
from the value that increases by 1 from the 
minimum value. Simple fitting error analysis of 
the flux density leads to very small errors (typi- 



2 



cally 1 mJy). We added flux errors to 10 % of the 
measured integrated flux density by assuming the 
typical accuracy of the absolute amp litude calibra- 
tion of the VLBA (|Cai et al. \\20m . 

The models for 2.3 GHz and 8.4 GHz are given 
in Tables 2 and 3, respectively. The fitted position 
of knots is denoted by a cross (+) and that of the 
core is denoted by an ellipse in the images. Com- 
ponent names in this paper are defined as follows: 
in the 2.3 GHz images, components A, B, C, D, 
and E, from downstream of the jet to the core Fs, 
and in the 8.4 GHz images, the most inner compo- 
nent E at 2.3 GHz is resolved into El, E2, E3, E4, 
and E5, from downstream of the jet to the core 
Fx. 

Knot identification at 2.3 GHz between Period 
1 and 2 is very clear, because no systematic motion 
of the knots can be found over our observation pe- 
riod. Time evolution of the position shift of each 
component appears to be random rather than sys- 
tematic, which implies that there is no evidence 
for proper motion in the 2.3 GHz images with 3 
mas resolution. In contrast, knot identification at 
8.4 GHz is not necessarily clear, because compo- 
nents El, E2, E3 and E4 appear to show expand- 
ing motion toward downstream of the jet. Thus, 
in this subsection, we first measure proper motion 
at 8.4 GHz in each Period separately. We then 
measure the combined proper motion between Pe- 
riod 1 and 2 in the next subsection, based on the 
results from the proper motion in each Period. In 
Tables 2 and 3, we have already referred to the 
identified names between Period 1 and 2, which 
are anticipated from the result presented in §3.2. 

To analyze the jet proper motion for each com- 
ponent, we carried out weighted least-squares fit- 
ting to the jet components detected at 8.4 GHz. 
Assuming that the core position is stationary, and 
that the jet moves along a straight line at an equal 
speed, we derived the fits to the R.A. direction (X) 
and Dec. direction (Y) separately. We then cal- 

+ u?, , and 



culated proper motion using u = 

position angle (PA) using (/)=tan~^(wx/uy), where 
Ux and Uy is the jet speed toward the X and Y 
directions, respectively. The results of our fits are 
shown in Table 4 and in Figures 5 and 6 for Period 
1 and Period 2, respectively. 

We estimated the proper motion of E3, E2, and 
El to be 0.21±0.01, 0.38±0.02, and 0.70±0.07 mas 



year ^, which corresponds to /3app= 0.30±0.02, 

0. 53±0.03, and 0.96±0.09, at 0.7, 1.5, and 2.5 mas 
(0.3, 0.7, and 1.1 pc) from the core during Period 

1, respectively. We estimated that of E4 and E3 
to be 0.23±0.03 and 0.31±0.02 mas year"!, which 
corresponds /3app= 0.32±0.04 and 0.43±0.03 at 0.7 
and 1.5 mas (0.3 and 0.7 pc) from the core dur- 
ing Period 2, respectively. The PA of the proper 
motion is different from that of the averaged jet 
direction at 8.4 GHz (57 degrees, see §3.4) by ap- 
proximately 30 degrees for E2 and 20 degrees for 
E3, implying possible curved trajectory of these 
knots. 

Dynamic range of our images at 8.4 GHz is not 
so high, typically 150. This would cause a domi- 
nant source of systematic position error in model 
fitting near the edge of the jet. In fact, the data 
points at the 2002.14 epoch and the 2004.81 epoch 
show a large deviation from the fitted line for the 
outer jet components (see Figures 5 and 6). 

3.2. Combined proper motion between 
Period 1 and Period 2 

In order to obtain the proper motion at 8.4 GHz 
over Period 1 and 2, a duration of more than 4 
years, we carried out knot identification between 
Periods 1 and 2. Figure 7 shows the combined 
proper motion of component E3. The motion Ux 
of component E3 in each Period is in good agree- 
ment with each other, which demonstrates that a 
constant speed is kept in the jet between 0.5 and 
1.6 mas (0.2 and 0.7 pc). However, the motion 
Uy of component E3 appears to be demonstrate 
acceleration, suggesting a curved trajectory. 

Although tracing of the jet trajectory is very 
interesting for the physics of jet propagation, we 
do not address it further in this paper because our 
observations did not cover over the whole period 
of the movement of E3. We simply fitted E3 over 
our whole observation period using a straight line 
as a first-order approximation, and obtained = 
0.245 ±0.002 mas year"!, Uy = 0.139 ±0.002 mas 
year-i, u = 0.282 ± 0.002 mas year"!, /3app = 
0.390 ± 0.003, and </) = 60.5 ± 0.6 degrees. 
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3.3. Calibration of the proper motion us- 
ing a model of orbital motion of the 
core 

The orbital motion of the core of 3C 66B with 
the orbital period of 1.05 ± 0.03 years has been 
shown by the pha se-referencing V LBI observations 



in Period 1 (S udou et al. 1 120031 ). This core shift 



could cause apparent position shifts of the jet com- 
ponents systematically, because we assumed that 
the core position is stationary in measuring the 
proper motion of the jet. This effect is expected 
to not be negligible if the position error of the 
knots is smaller than the major axis of the orbital 
motion, which was estimated to be 0.080 mas at 
8.4 GHz. 

Assuming that the core position in our images 
moves in the same way as the core orbital motion, 
we are able to calibrate for the effect of the core 
motion and to expect to obtain a more accurate 
measurement of the proper motion of the jet. The 
core orbital motion is modeled on the basis of the 
presence of a binary black hole with a circular or- 
bit in the core of 3C 66B, a method which was 



applied in lSudou et al. I (|2003l ). The obtained cal- 
ibration table is presented in Table 5. 

We show the calibrated proper motion of the jet 
in Figures 8 and 9 and Table 6 in the same manner 
as Figures 5 and 6 and Table 4. Significant differ- 
ences between the calibrated and non-calibrated 
proper motion could not be found. The fitting re- 
sult of the most inner component, E3, gave a much 
better xt value than that of the non-calibrated 
one in both the X and Y directions. Calibration 
toward the X direction at the 2002.45 epoch will 
work very effectively, because the core position 
at this epoch is isolated from those at the other 
epochs indicated in Table 5. Indeed, the decrease 
in the xt value of E3 is mainly caused by the po- 
sition shift at the 2002.45 epoch. The other com- 
ponents exhibited xt values that were unchanged 
or worse in either or both directions. It is nat- 
ural that the calibration of the outer jet compo- 
nent. El, with a position error larger than the core 
orbital diameter (0.08 mas) would be inefficient. 
Calibration for components E3 and E4 in Period 
2 also appear to be inefficient, even though their 
position errors are smaller than the orbital radius. 
Originally, their non-calibrated fit of showed 
the xt value much larger than unity (58.4 and 



45.0, respectively), indicating their position shift 
contains some unexpected error. This is probably 
why calibration is not effective for these compo- 
nents. 

3.4. Structure of the jet 

The images at 2.3 GHz show the straight and 
smooth jet structure with several knots and the 
unresolved core in 3C 66B at scales within a dis- 
tance of 20 mas (9 pc) from the nucleus (Figures 
1 and 2). The jet extends toward the north-east 
direction with the PA of 55 ± 4 degrees, which 
is estimated from the distribution of the knots. 
Weak wiggling patterns can be found, in particular 
at epoch 2001.20, 2001.48, 2002.14, and 2002.45. 
These parsec-scale characteristics agree with re- 
sults in the VLBI ima ges previously obtained at 



1.6 GHz and 5.0 GHz ( Fraix-Burnet et al. 1 11997 



Xu et al. |[2000t iGiovannini et al. II2OOII) 



The images at 8.4 GHz show the sub-parsec 
scale jet structure at scales within a distance of 
3 mas (1.3 pc) from the nucleus (Figures 3 and 
4). At this frequency, the jet becomes relatively 
weaker and smoother than at 2.3 GHz. The aver- 
aged PA of the jet is estimated to be 57 ± 4 de- 
grees. The difference of the PA between 2.3 GHz 
and 8.4 GHz is only a few degree, indicating a 
well-aligned jet structure from sub-parsec to par- 
sec scales. A weak bending jet pattern toward the 
east direction at a distance of 2 mas from the core 
is seen in 8.4 GHz images, in particular at epoch 
2001.20, 2001.48, 2002.14, and 2005.69 (Figures 3 
and 4). 

We were unable to find the co unterjet that was 
detect ed in the 5.0 GHz image bv lGiovannini et al. 
(|200l[ ) . Assuming that the reported jet/ counterjet 
ratio of 10 can be applied to our images, we ex- 
pect the counterjet intensity to be 0.5 mJy which 
is smaller than the 3-(T noise limit of our images 
at both frequencies. 

3.5. Flux variation of the jet and the core 

Time evolution of the flux density of the jet 
components are shown in Figure 10. Note that 
the sudden flux dip at the 2005.35 epoch at both 
frequencies may be affected by the lower resolu- 
tion due to the lack of longer baseline data related 
to the St. Croix station. All the components seem 
to stay at an unchanged flux level within measure- 
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ment errors. By assuming the jet flux densities are 
constant, we can obtain an absolute amplitude er- 
ror of 15 %. This is consistent with the typical 
error that we assumed in §2. 

We also found that the flux density of the core 
increased by about 50 % at 8.4 GHz (see Fig- 
ure 11). This flux curve of the core is uncorre- 
lated with those in the je t comp onents. Recently 

(|2010l ) revealed the 93- 
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day variation of the mm- wave flux density of the 
core region, and interpreted it as periodic Doppler 
modulation due to the orbital motion of a sub- 
parsec scale binary black hole. This variation is 
expected to be hard to be detected at cm-waves, 
because Syncrotron self- absorption would be dom- 
inated in the nucleus at lower radio frequencies 
and dilute the variation amplitude. Unfortunately 
the time sampling of our observations is too sparse 
and uneven (10 to 100 days), this short-term vari- 
ation cannot be discussed in this paper. 

3.6. Structural change of the core 

Figure 11 shows the changes of the integrated 
flux density, major axis, and PA of the fitted pa- 
rameters of the core. At 8.4 GHz, we found that 
the major axis increased with time from 0.28 to 
0.38 mas in Period 1, and that after reaching its 
largest value of 0.42 mas, decreased to 0.32 mas 
at the last epoch in Period 2. The PA changed 
from 80 to 50 degrees in Period 1. A possible in- 
crease of the flux density and major axis can be 
also seen at 2.3 GHz. This extension of the major 
axis and flux enhancement of the core can be ex- 
plained qualitatively by a new component emerg- 
ing that propagates in the core region. 

The axial ratio of the core at both frequencies 
equals to zero at almost all epochs, indicating that 
the core cannot be resolved toward perpendicular 
to the jet axis on the sub-parsec scale. 

4. Discussion 

We plotted the apparent speed of the jet as 
a function of the distance from the core in Fig- 
ure 12a. 3C 66B shows that the apparent jet 
speed increases with a distance from the core 
within 3 mas (1.3 pc). This result agrees with 
the fact that the sub-parsec scale counterjet was 
found within a distance of 2 mas from the core 
( Giovannini et al. I 120011 ). Since the measured 



proper motion in 3C 66B is considered to reflect 
the bulk speed of the jet within ob servational er- 
ror (jlguchi. Okuda. fc Sudou1l2010[ ). the apparent 
jet speed can be written as 



/3app — 



13 sin i 
1 — /3 cos i ' 



(1) 



where /3 is the jet outflow speed (=v/c) and i is the 



viewing angle of the jet (|Ginzburg fc Svrovatskii 
19691) . This equation indicates that increases in 



the apparent speed are caused by either jet ac- 
celeration or a change in the viewing angle. Fig- 
ure 12b shows the constraints of the jet speed 
and t he viewing angle accord i ng to Equation 
(1). Ilguchi. Okuda. fc Sudou! (|2010h revealed 
that /3 = 0.77 and i = 5.3 degrees using the mm- 
wave periodic variation together with the proper 
motion of E3 in Period 1. 

Assuming that the jet accelerates with a con- 
stant viewing angle of 5.3 degrees, the jet speed 
accelerates from 0.77 c to 0.92 c, correspond- 
ing to the bulk Lorentz factor F from 1.60 to 
2.57. Since the relativistic energy can be writ- 
ten as E = Tmv^, the kinetic energy gain by 
this acceleration corresponds to 1 GeV for normal 
plasma and 0.6 MeV for pair plasma. This energy 
gain may relate to a possible 7-ray emission from 
3C66B. Interestingly a very high energy (VHE) 7- 
ray emission above 150 GeV was detected in the 
region that includes both 3C 6 6B and a blazar 3C 
66A by the MAGIG tel e scope (|Aliu et al. II2009I ). 



Tavecchio fc Ghisellini ( 20091) discussed a possi- 
ble VHE 7-ray mechanism from 3C 66B in an anal- 
ogy of blazar emission models. If indeed 3C 66B 
is a VHE source, the VHE 7-ray emission might 
be related to the sub-parsec scale acceleration of 
the jet. However, we cannot discuss any more the 
detailed mechanism of the 7-ray emission from 3C 
66B, because physical conditions in the nuclear re- 
gion are unclear from very limited SED data, e.g., 
a thermal protons density for proton-proton in- 
teraction or high-energy electrons density for the 
inverse-Compton mechanism. 

Assuming the other interpretation, that is, a 
change in the viewing angle with a constant jet 
velocity of 0.77 c, the angle varies from 5.3 to 20.8 
degrees. This interpretation is supported from the 
presence of the weak apparent bend in the jet 
at a distance of 2 mas from the core described 
in §3.4. However, there are some sources which 
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show an apparent increase in the jet speed with 
a distance from pc to 100-pc scales without an 
bends in the jet, e.g., C en A (jTingav et al. Ill99 



Hardcastle et al. II2003I) . This may indicate that 



the relativistic jet model expressed in equation (1) 
seems to be too simple to explain general proper- 
ties of the apparent jet acceleration. We may need 
to consider more sophisticated models, for exam- 
ple, changes of emissivity profile with the distance 
from the core in the stratified (two-zone) jet with 
the spine and layer structure. In this case, it is 
expected that the spine jet is highly relativistic 
and that the layer jet is mildly relativistic. If the 
two-zone jet model is correct, the speed of each 
component would be discontinuously changed at 
the boundary between the region where emission 
from the layer jet is dominated and the region 
where emission from the spine jet is dominated. 
In case of 3C66B, the boundary is expected to be 
the distance of around 2 mas from the core. On 
the other hand, if the single-zone jet model which 
we assumed in equation (1) is correct, the speed 
of each component would be increased gradually 
with the distance from the core. For this pur- 
pose, a longer-span observation of several years is 
needed to trace the trajectory of the individual 
component in the parsec-scale jet of 3C 66B. 

5. Conclusion 

We have presented VLBI images at 2.3 GHz and 
8.4 GHz at 10 epochs over 4 years. Subluminal 
motion of the jet ranging from 0.30 c to 0.96 c was 
detected within a distance of 3 mas (1.3 pc) from 
the core from the 8.4 GHz images. We demon- 
strated the calibration of the core position using 
the result of phase-referencing VLBI observations, 
which could improve the accuracy of measuring 
subluminal proper motion of a sub-parsec scale jet 
in nearby radio galaxies. Using a model of the core 
motion due to an orbital motion of a binary black 
hole found in 3C 66B, we calibrated the proper 
motion of the jet, and found no significant dif- 
ference between the calibrated and non-calibrated 
proper motions. Our observational results indicate 
that the proper motion at sub-parsec scale is likely 
to be increased by changes in the viewing angle, 
for instance, the angle varies from 5 to 21 degrees 
in case of a constant jet speed of 0.77 c. However, 
the possibility of acceleration of the jet outflow 
speed still remains. Another interpretation is also 



suggested that changes of emissivity profile in the 
two-zone jet cause the changes of the proper mo- 
tion. To reveal the kinematics and geometry of the 
jet in 3C 66B, further monitoring is required over 
a longer period of several years to trace the indi- 
vidual jet components from sub-parsec to parsec 
scales. 
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Table 1 

Observations and image parameters 



iL/pocn 


rreq. 


Array 


Int. lime 


Beam Size 


Beam PA 


rms noise 








[hrj 


[mas X mas] 


[dcg] 


[m Jy /beam] 








(A \ 




(b) 


{<) 


2001.20 


2.3 


BFHKLMNOPS 


0.7 


4.0 X 2.5 


-1.7 


0.5 


2001.48 


2.3 


BFHKLMNOPS 


0.3 


4.0 X 2.6 


9.4 


0.6 


2001.86 


2.3 


BFHKLMNOPS 


0.4 


3.7 X 2.8 


-9.7 


0.6 


2002.11 


2.3 


BFHKLMNOPS 


0.4 


4.0 X 2.7 


-7.7 


0.4 


2002.14 


2.3 


BFHKLMNOPS 


0.4 


4.0 X 2.7 


-5.2 


0.6 


2002.45 


2.3 


BFHKLMNOPS 


0.4 


4.0 X 2.7 


-5.0 


0.5 


2004.80 


2.3 


BFHKLMNOPS 


0.4 


3.9 X 2.4 


6.8 


0.7 


2005.05 


2.3 


BFHKLMNOPS 


0.5 


3.8 X 2.4 


9.2 


0.7 


2005.35 


2.3 


BFHKLMNOPS^' 


0.4 


6.4 X 3.5 


-6.5 


0.7 


2005.54 


2.3 


BFHKLMNOPS 


0.5 


3.8 X 2.3 


3.2 


0.7 


2001.20 


8.4 


BFHKLMNOPS 


0.6 


1.0 X 0.6 


1.4 


0.8 


2001.48 


8.4* 


BFHKLMNOPS 


0.3 








2001.86 


8.4 


BFHKLMNOPS 


0.4 


1.1 X 0.7 


-8.6 


0.8 


2002.11 


8.4 


BFHKLMNOPS 


0.4 


1.0 X 0.7 


-7.0 


0.9 


2002.14 


8.4 


BFHKLMNOPS 


0.4 


1.0 X 0.7 


-8.5 


0.6 


2002.45 


8.4 


BFHKLMNOPS 


0.4 


1.1 X 0.7 


-3.3 


0.7 


2004.80 


8.4 


BFHKLMNOPS 


0.5 


1.1 X 0.6 


4.8 


1.2 


2005.05 


8.4 


BFHKLMNOPS 


0.4 


1.0 X 0.6 


-2.0 


0.7 


2005.35 


8.4 


BFHKLMNOPS^^ 


0.4 


1.4 X 0.6 


-21.0 


1.0 


2005.54 


8.4 


BFHKLMNOPS 


0.5 


0.9 X 0.6 


-3.0 


0.8 


2005.69 


8.4 


BFHKLMNOPSE 


1.7 


0.8 X 0.6 


18.4 


0.4 



Note. — (1) Observing Epoch, (2) Observing Frequency, (3) Array configuration; 
B=Brcwstcr, F=Fort Davis, H=Hancock, K=Kitt Peak, L=Los Alamos, M=Maiina Kca, 
N=North Liberty, 0=Owens Valley, P=Pie Town, S=St. Croix, E=Efl[elsberg, the diame- 
ter of the antennas is 100 m for Effelsberg and 25 m for the VLBA antennas, (4) Integration 
time, (5) Major axis and minor axis of the CLEAN beam, (6) Position angle of the CLEAN 
beam, (7) rms noise level of image. 

* We removed this data at 8.4 GHz from imaging analysis because of low quality of the data 
(see §2). 

^ The data of the St. Croix station was completely removed because of instrumental problems. 
This makes the beam size about 1.5 times larger than that of the full VLBA configuration. 
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Table 2 

Positions and flux densities of the knots at 2.3 GHz 



Epoch 


Comp. 


F 


Ax 


Ay 


Major axis 


Axial 


PA 






r T 1 

[Jy] 


[masj 


[masJ 


[masJ 


ratio 


[deg] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


JrtJiiUU. 1 
















9001 90 


s 


1 ni 9 -1- n 01 09 


n 00 -1- 00 


n 00 -1- 00 
u.uu m u.uu 


1 7fi -1- 09 
1. / zlz U.UZ 


00 -1- 00 
u.uu ZlZ U.UU 


oy.o ziz u.o 


iXi/— ^-^ ) 




01 1 -1- 001 A. 

u.u±±u m u.uuj-^ 


Q QQ -1- 0"? 


1 CO -1- A A/l 

±.oz in u.yi 










D 


n 0048 ± 0009 


5 24 ± 08 


3 84 ± 07 










n 


n nO"^? + OOOQ 


O.WJ —1— U.-L-L 


t; CA -1- A 1 1 










J3 


n 0048 + OOOQ 


1 fiK _(_ A nq 


q 71 -1- A 1 A 

ij • i L _l_ U. J-U 








9001 Ak 


-T s 


n 1 079 -1- 01 Of^ 

U.IU/Z HZ U.UIUO 


n 00 -1- 00 

u.uu ZIZ u.uu 


n 00 -1- 00 
u.uu ZIZ u.uu 


1 9Q _|_ n 09 

L.ZiO ZIZ U.UZ 


n 00 -1- 00 

u.uu ZIZ u.uu 


(:;7 9 _|_ 9 n 

U 1 . Z ZIZ Zd.xj 


IXi/— ■'■•o ) 




n 01 1 R -1- 001 ^ 


9 _l_ n OA 
z.yo ZIZ u.U'4: 


9 OR -1- 04 

Z.UO ZIZ U.U'4: 










D 


0044 ± 0011 


5.72 ± 0.11 


3.44 ± 0.11 










n 


n nnQs + n noi 


q (-)» + 1 4 


G QO _|_ A 1 4 










B 


n 0047 + 001 1 


11 07 -I- n 10 


q 21 +0 12 












n 00*^4 + 001 


i« ni -1- n 14 


1 Q c;a _|_ A 1 








9nni sfi 


^ s 


n 1 097 -1- 01 0*^ 

U.lUZff ZIZ U.UIUO 


fl 00 -1- 00 
u.uu ZIZ u.uu 


n OzL -1- OA 

U.U4: ZIZ U.U^ 


1 _|_ n 09 

1.00 ZIZ u.uz 


n 00 -1- 00 
u.uu ZIZ u.uu 


00.0 ZIZ l.fl 


9 "1 


11/ 


n ni '^'^ -1- 001 

U.UIOO ZIZ U.UUlO 


9 00 _|_ n fi9 

Zi.OO ZIZ U.UZi 


1 Qo _|_ n 09 
±.yu ZIZ u.uz 










D 


0057 ± 0008 


5 38 ± 05 


4 04 ± 05 










V, 


n 0057 + 0008 


q Q7 _(_ A QG 


« GG _(_ A AG 

U.UU —1— u.uu 










B 


n 005Q + 0008 


XZ(.4:U —1— U.UU 


1 A 40 -1- A A« 

xu.^ic/ —1— u.uu 










Pc. 
-T S 


n 1 91 1 _|_ n 01 91 

U.J-Z(_L_L ZIZ U.UIZI 


n 00 -1- 00 

u.uu ZIZ u.uu 


n 00 -1- 00 
U.UU ZIZ u.uu 


1 ztR -1- 09 

L.'^O ZIZ U.UZ 


00 -1- 00 
u.uu ZIZ u.uu 


Urt.O ZIZ ±.J- 


\Xi^— ^-^ ) 




U.UIOI ZIZ U.UUlU 


9 r^7 _|_ [\ no 
z.u i ZIZ u.uo 


1 _|_ A MO 

i.uo ZIZ u.uo 












n 00Q9 -1- 001 '\ 

U.UUc/Z ZIZ U.UUlO 


'i.oo ZIZ u.uu 


00 _|_ n fif^ 
0.00 ZIZ u.uu 












00(^/1 -U 001 1 
U.UUD4 zt U.UUii 


Q 77 _|_ 1 7 
o.i i zt U.i / 


CK QO -U 00 

D.ou zt u.uy 










B 


0.0065 ± 0.0011 


12.35 ± 0.08 


10.04 ± 0.09 










A 


0.0036 ± 0.0010 


14.98 ± 0.17 


13.53 ± 0.18 








2002.14 


Fs 


0.1229 ± 0.0123 


0.00 ± 0.00 


0.00 ± 0.00 


1.58 ± 0.02 


0.19 ± 0.19 


51.3 ± 0.8 


{Xl= 2.1 ) 


E 


0.0132 ± 0.0017 


2.59 ± 0.03 


1.66 ± 0.03 










D 


0.0076 ± 0.0013 


4.39 ± 0.06 


3.35 ± 0.06 










C 


0.0051 ± 0.0011 


7.00 ± 0.17 


5.56 ± 0.15 










B 


0.0062 ± 0.0012 


12.21 ± 0.08 


9.98 ± 0.09 










A 


0.0028 ± 0.0010 


15.75 ± 0.17 


14.73 ± 0.19 









Table 2 — Continued 



Epoch 


Comp. 


F 


Ax 


Ay 


Major axis 


Axial 


PA 






[Jy] 


[mas] 


[mas] 


[mas] 


ratio 


[deg] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




2002.45 


Fs 


0.1137 ± 0.0114 


0.00 ± 0.00 


0.00 ± 0.00 


1.64 ± 0.02 


0.00 ± 0.00 


59.6 ± 1.0 


(y'^= 2.1 ) 


E 


0.0123 ± 0.0015 


2.85 ± 0.03 


2.06 ± 0.04 










D 


0.0056 ± 0.0011 


4.93 ± 0.08 


2.99 ± 0.08 










c 


0.0068 ± 0.0012 


8.40 ± 0.07 


6.12 ± 0.07 










B 


0.0063 ± 0.0011 


11.72 ± 0.09 


9.14 ± 0.09 










A 


0.0031 ± 0.0010 


15.08 ± 0.14 


13.29 ± 0.16 








Period 2 
















2004.80 




0.1108 ± 0.0111 


0.00 ± 0.00 


0.00 ± 0.00 


1.69 ± 0.01 


0.00 ± 0.00 


64.4 ± 1.1 


(■y^= 2.1 ) 


E 


0.0121 ± 0.0015 


2.87 ± 0.04 


2.12 ± 0.04 










D 


0.0047 ± 0.0010 


4.99 ± 0.10 


2.86 ± 0.10 










C 


0.0059 ± 0.0011 


8.44 ± 0.08 


6.02 ± 0.08 










B 


0.0045 ± 0.0010 


11.84 ± 0.13 


9.42 ± 0.13 








2005.05 


Fc; 


0.1011 ± 0.0101 


0.00 ± 0.00 


0.00 ± 0.00 


1.70 ± 0.02 


0.00 ± 0.00 


58.7 ± 0.9 


(y2= 2.0 ) 


E 


0.0122 ± 0.0014 


2.71 ± 0.03 


1.49 ± 0.03 










D 


0.0044 ± 0.0008 


4.51 ± 0.08 


3.35 ± 0.08 










c 


0046 + OOOQ 


7 49 + 08 


5 21 + 08 










B 


0.0053 ± 0.0009 


11.41 ± 0.09 


8.76 ± 0.08 








2005.35 


Fs 


0.1157 ± 0.0116 


0.00 ± 0.00 


0.00 ± 0.00 


2.36 ± 0.04 


0.00 ± 0.00 


50.9 ± 1.3 


{Xl= 2.0 ) 


E 


0.0067 ± 0.0013 


2.72 ± 0.09 


1.99 ± 0.09 










D 


0.0086 ± 0.0014 


4.46 ± 0.11 


3.00 ± 0.12 










C 


0.0057 ± 0.0012 


8.59 ± 0.10 


6.09 ± 0.11 










B 


0.0074 ± 0.0013 


12.68 ± 0.12 


9.42 ± 0.13 








2005.54 


Fs 


0.1271 ± 0.0127 


0.00 ± 0.00 


0.00 ± 0.00 


1.69 ± 0.02 


0.00 ± 0.00 


53.9 ± 0.9 


{Xi= 2.3 ) 


E 


0.0162 ± 0.0019 


2.92 ± 0.03 


1.63 ± 0.03 










D 


0.0059 ± 0.0011 


5.10 ± 0.08 


3.10 ± 0.08 










C 


0.0060 ± 0.0011 


8.71 ± 0.09 


5.95 ± 0.09 










B 


0.0060 ± 0.0011 


12.00 ± 0.10 


9.03 ± 0.10 









Note. — (1) Observing Epoch, (2) Component names defined as follows; components A, B, C, D, and E, from 
downstream of the jet to the core Fg, (3) Flux density, (4) R.A. offset, (5) Dec. offset, (6) Major axis of Gaussian 
components, (7) Axial ratio of Gaussian components, (8) Position angle of Gaussian components. 



Table 3 

Positions and flux densities of the knots at 8.4 GHz 



Epoch 


Gomp. 


F 


Ax 


Ay 


Major axis 


Axial 


PA 






[Jy] 


[mas] 


[mas] 


[mas] 


ratio 


[deg] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


PprioH 1 
















9001 90 


^ X 


n OSS + nriQ 




A nn _|_ n AA 
u.uu zn u.uu 


A OQ -1- A ni 

u.^o zn U.UJ. 


n no + n 00 

u.uu zn u.uu 


7Q n + 1 fi 

( o.u zn J. .u 


{Xv— -L-U ) 




n ni c _(_ ri nno 


n cri _(_ ri ni 

U.tJU —1— U.UJ- 


A 40 _|_ n Ai 










E2 


n r)r)7 + n nni 


1 1 + 09 

J..±U —1— u.u^ 


n QG _|_ n no 

U.OU —1— U.UZi 










El 


0.005 ± 0.001 


1.85 ± 0.05 


1.03 ± 0.05 








^uux .ou 


X 




n no + 00 
u.uu zn u.uu 


00 + 00 
u.uu zn u.uu 


n qi -(- n 01 

u.OJ. zn U.UJ. 


n no + 00 

u.uu zn u.uu 


CO Q _(_ 1 f; 
00.0 zn J..U 


{Xv— u.y ) 




n 01 8 + 009 


n etc _|_ ri ni 


A 44 _|_ n ni 

U.t:t: —I— U.U J- 










E2 




1 QO _(_ A no 


A oq _(_ n AO 










El 


0.003 ± 0.001 


2.13 zb 0.07 


1.22 zb 0.07 








9009 1 1 

ZiWW^ • -L J. 


Fv 
X 


n 1 HQ _(_ n ni 1 


n nn _(_ n nn 
U.UU —1— u.uu 


A nn _(_ n AA 
U.UU —1— u.uu 


n ^7 + n 01 

\J.K> I —1— U.U-L 


n 1 q _(_ n ns 

\J.±^ _l_ u.uo 


7 _(_ Q n 

uo. I —1— o.u 






n ni 7 -1- n nno 


n 71 -\- ni 

U . i -L —1— U.U -L 


A 41 + n 01 

U.t:-L —I— U.UJ- 










F,9 


n nrift 4- 001 


1 44 09 

ZIZ U.UZi 


n 04 _L n no 

U.O^ ZIZ U.UZi 










El 


0.003 ± 0.001 


2.43 zb 0.08 


1.45 zb 0.08 










Fv 
^ X 


n 1 1 q -1- n ni o 


n nn -1- n nn 
U.UU —1— u.uu 


n nn -1- n nn 
U.UU —1— u.uu 


n Q4 + n 01 

U.t>*± —1— U.U-L 


08 + 08 

U.UO —1— u.uo 


C _|_ Q 1 
UO.U —1— O.-L 






n ni q _(_ n nno 


n 70 -1- n ni 
yj , I Zj _i_ U.U X 


A 40 + n AI 










E2 


n nnc _(_ n nni 

U.UUO —1— u.uux 


1 Q7 _|_ A no 


f) 00 1 n no 










El 


n nno _l n nni 

U.UUO —1— u.uux 


on _(_ n nq 


n nc. _L n nn 

U.t/U —1— U.Ut? 








2002.45 




127 ± 013 


00 ± 00 


00 ± 00 


38 zb 02 


00 zb 00 


48.7 zb 4.1 






n ni o _L n nno 

U.U-LO —1— U.UUZi 


n 71 _L no 


f) AO _L n no 










E2 


n nin -I- n nni 

U.UJ-U —1— U.UU-L 


1 «i -I- A no 


A oc 1 n no 

U.OU —1— U.UZj 










F1 


f\f\A -\- n 001 


9 7c _l_ n nn 


1 c^n -1- HQ 

i.ou ZIZ U.uy 








Period 2 
















2004.80 


Fx 


0.095 ± 0.010 


0.00 zb 0.00 


0.00 zb 0.00 


0.39 zb 0.02 


0.00 zb 0.00 


51.6 zb 3.3 


(X^= 1-9 ) 


E4 


0.018 ± 0.002 


0.78 zb 0.02 


0.49 zb 0.02 










E3 


0.007 ± 0.001 


1.58 ± 0.02 


0.99 ± 0.02 










E2 


0.005 ± 0.001 


3.04 zb 0.09 


1.91 zb 0.08 








2005.05 


Fx 


0.092 ± 0.009 


0.00 zb 0.00 


0.00 zb 0.00 


0.42 zb 0.01 


0.00 zb 0.00 


54.0 zb 3.3 


{Xl= 1-9 ) 


E4 


0.018 ± 0.002 


0.68 zb 0.02 


0.51 ± 0.02 










E3 


0.010 ± 0.001 


1.39 zb 0.01 


0.85 zb 0.01 









Table 3 — Continued 



Epoch 
(1) 



2005.35 

{xi= 2.0 



Comp. 

(2) 



F 

[Jy] 

(3) 



Ax 
[mas] 
(4) 



Ay 
[mas] 
(5) 



Major axis 
[mas] 
(6) 



Axial 
ratio 
(7) 



PA 
[deg] 
(8) 



Fx 
E4 
E3 



0.091 ± 0.009 
0.012 ± 0.002 
0.007 ± 0.001 



0.00 ± 0.00 
0.71 ± 0.03 
1.57 ± 0.03 



0.00 ± 0.00 
0.59 ± 0.03 
1.07 ± 0.03 



0.33 ± 0.02 0.00 ± 0.00 46.0 ± 6.8 



2005.54 

{xl= 1.9 ; 



Fx 
E4 

E3 



0.120 ± 0.012 
0.020 ± 0.002 
0.010 ± 0.001 



0.00 ± 0.00 
0.69 ± 0.02 
1.47 ± 0.02 



0.00 ± 0.00 
0.48 ± 0.02 
0.94 ± 0.02 



0.40 ± 



01 0.13 ±0. 



07 57.8 ± 



2005.69 

{xi= 1-2 : 



Fx 
E5 
E4 
E3 



0.111 ± 0.011 
0.020 ± 0.002 
0.010 ± 0.001 
0.006 ± 0.001 



0.00 ± 0.00 
0.53 ± 0.02 
1.00 ± 0.02 
1.66 ± 0.01 



0.00 ± 0.00 
0.31 ± 0.02 
0.70 ± 0.02 
1.02 ± 0.01 



0.32 ± 



01 0.00 ± 0. 



00 56.2 ± 



Note. — Same as Tabic 2, but at 8.4 GHz. Component names defined as follows; components El, E2, E3, 
E4 and E5 from downstream of the jet to the core Fx- 



Table 4 

Proper motion of the jet in 3C66B 



Comp. 




Uy 


u 


^app 


PA 




[mas / yr] 


[mas/vr] 


[mas / vr] 




[dcg] 


(1) 


(2)' 




(4)" 


(•"') 


(0) 


Period 1 












El 


0.60 ± 0.07 


0.35 ± 0.08 


0.70 ± 0.07 


0.96 ± 0.09 


60 ± 8 


E2 


0.38 ± 0.02 


-0.02 ± 0.02 


0.38 ± 0.02 


0.53 ± 0.03 


94 ± 3 


E3 


0.21 ± 0.01 


0.01 ± 0.01 


0.21 ± 0.01 


0.30 ± 0.02 


88 ± 3 


Period 2 












E3 


0.27 ± 0.02 


0.16 ± 0.02 


0.31 ± 0.02 


0.43 ± 0.03 


59 ± 5 


E4 


0.17 ± 0.03 


0.16 ± 0.03 


0.23 ± 0.03 


0.32 ± 0.04 


46 ± 10 



Note. — (1) Component names defined as follows; components El, E2, E3, and 
E4 from downstream of the jet to the core Fx, (2) Jet proper motion toward R.A. 

direction, (3) Jet proper motion toward Dec. direction, (4) Jet proper motion, (5) 
Apparent jet velocity /3app = I'app/c, (6) Position angle of the proper motion 
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Table 5 

Calibration table of the core position shift using a model of a binary black hole 



Epoch 


Ax 


Ay 


[year] 


[mas] 


[mas] 


(1) 


(2) 


(3) 


Period 1 






2001.20 


-0.010 


0.010 


2001.86 


-0.024 


-0.009 


2002.11 


-0.040 


0.005 


2002.14 


-0.035 


0.006 


2002.45 


0.038 


0.006 


Period 2 






2004.80 


0.027 


-0.008 


2005.05 


-0.034 


-0.008 


2005.35 


-0.026 


0.008 


2005.54 


0.025 


0.009 


2005.69 


0.045 


0.002 



Note. — (l)Observation 
epoch, (2) R. A. offset from 
the center of the core orbit, (3) 
Dec. offset from the center of 
the core orbit. These param- 
eters are calculated based on 



Table 1 of ISudou et al. I (120031) 
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Table 6 

Proper motion of the jet in 3C66B calibrated by a model of core motion 



Comp. 




Uy 


u 


/3app 


PA 




[mas/yr] 


[mas/yr] 


[mas/yr] 




[deg] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Period 1 












El 


0.61 ± 0.07 


0.35 ± 0.08 


0.70 ± 0.07 


0.96 ± 0.09 


61 ± 8 


E2 


0.40 ± 0.02 


-0.02 ± 0.02 


0.40 ± 0.02 


0.55 ± 0.03 


95 ± 3 


E3 


0.20 ± 0.01 


0.01 ± 0.01 


0.20 ± 0.01 


0.28 ± 0.02 


88 ± 3 


Period 2 












E3 


0.34 ± 0.02 


0.17 ± 0.02 


0.39 ± 0.02 


0.53 ± 0.03 


63 ± 4 


E4 


0.20 ± 0.03 


0.18 ± 0.03 


0.27 ± 0.03 


0.37 ± 0.04 


49 ± 8 



Note. — (1) Component names defined as follows; components El, E2, E3, 
and E4 from downstream of the jet to the core Fx, (2) Jet proper motion toward 
R.A. direction, (3) Jet proper motion toward Dec. direction, (4) Jet proper 
motion, (5) Apparent jet velocity /3app = Vapp/c, (6) Position angle of the proper 
motion. 
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Table 7 

values of the linear fitting of proper motions before and after applying the 

calibration for the core motion. 



Comp. X-direction Y-direction 





N 


C 


N 


C 


(1) 


(2) 


(3) 


(4) 


(5) 


Period 1 










El 


2.1 


4.0 


3.5 


5.0 


E2 


6.4 


16.2 


1.2 


0.7 


E;] 


4.1 


0.1 


3.() 


3.2 


Period 2 










E3 


58.4 


91.0 


30.2 


31.2 


E4 


45.0 


62.0 


16.7 


15.0 



Note. — (1) Component name de- 
fined in Table 4, (2) xl value for R.A. di- 
rection without the correction of the core 
orbital motion, (3) xl value for R.A. di- 
rection with the correction of the core 
orbital motion, (4) xt value for Dec. di- 
rection without the correction of the core 
orbital motion, (4) xt value for Dec. di- 
rection with the correction of the core or- 
bital motion. Only the most inner com- 
ponent, E3, exhibited a better xl value 
in both the R.A. and Dec. directions, 
the others exhibited unchanged or worse 
xl values in either or both directions. 
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Right Ascension (mas) 



Fig. 1.— Six 2.3 GHz images of 3C 66B during 2001.20 to 2002.45 (Period 1). Contours are drawn at 
3-cr X (1, V2,2,2\/2, ...). Each value of a is 0.5, 0.6, 0.6, 0.4, 0.6 and 0.5 mJy beam"! from (a) to (f), 
respectively. The synthesized beam is shown at the bottom-left of the image, the accurate size of which is 
also shown in Table 1. The cross marks indicate the fitted position of knots using a delta functionin. The 
Gaussian component of the core is also shown. 
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Fig. 2.— Four 2.3 GHz images of 3C66B during 2004.80 to 2005.54 (Period 2). Contours are drawn at 
3-tT X (1, v^, 2, 2-\/2, ...)• Each value of a is 0.7 mJy beam"^ for these epochs. The synthesized beam is 
shown at the bottom- left of the image, the accurate size of which is also shown in Table 1. The cross marks 
indicate the fitted position of knots using a delta functionin. The Gaussian component of the core is also 
shown. 
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Fig. 3.— Five 8.4 GHz images of 3C66B during 2001.20 to 2002.45 (Period 1). Contours are drawn at 3- 
a X (1, v^, 2, 2\/2, ...). Each value of a is 0.8, 0.8, 0.9, 0.6 and 0.7 mJy beam^^ from (a) to (e), respectively. 
The synthesized beam is shown at the bottom-left of the image, the accurate size of which is also shown 
in Table 1. The cross marks indicate the fitted position of knots using delta functionin. The Gaussian 
component of the core is also shown. 
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Fig. 4.— Five 8.4 GHz images of 3C66B during 2004.80 to 2005.69 (Period 2). Contours arc drawn at 3- 
cr X (1, v^, 2, 2-\/2, ...). Each value of u is 1.2, 0.7, 1.0, 0.8 and 0.4 mJy beam^^ from (a) to (e), respectively. 
The synthesized beam is shown at the bottom-left of the image, the accurate size of which is also shown 
in Table 1. The cross marks indicate the fitted position of knots using a delta functionin. The Gaussian 
component of the core is also shown. 
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Fig. 5. Linear fitting of tfic relative position shift of the knots El, E2, and E3 at 8.4 GHz in Period 1. 
(a) time evolution and (b) spatial distribution. A filled lozenge means the R.A. direction (X), and an open 
lozenge means the Dec. direction (Y). Dotted lines indicate the best fit. 
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Fig. 7. Combined proper motion of the component E3 between Period 1 and Period 2. Dotted lines 
indicate the fitting result in each Period, which was also shown in Figures 5 and 6, and the dashed line 
indicates the linearly fitted line of combined data of E3 between Period 1 and Period 2. 
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Fig. 8. Same as Figure 5, but the core position is ealibrated by a model of circular orbital motion in a 
putative binary black hole in 3C 66B. (a) time evolution and (b) spatial distribution. A filled lozenge means 
the R.A. direction (X), and an open lozenge means the Dec. direction (Y). Dotted lines indicate the best fit. 
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Fig. 9. — Same as Figure 6, but the core position is calibrated by a model of circular orbital motion in a 
putative binary black hole in 3C 66B. (a) time evolution and (b) spatial distribution. A filled lozenge means 
the R.A. direction (X), and an open lozenge means the Dec. direction (Y). Dotted lines indicate the best fit. 
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Fig. 10. — Flux variation of the jet components, (a) at 2.3 GHz, and (b) at 8.4 GHz. The components A at 

2.3 GHz and E5 at 8.4 GHz are excluded in this figure, because of insufficient number of data points. No 
significant systematic variation cannot be found at both frequencies. Note that the sudden flux dip at the 
2005.35 epoch at both frequencies may be affected by the lower-resolution due to the lack of longer baseline 
data related to the St. Croix station. 
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Fig. 11. — Time evolution of the fitted core parameter (a) at 2.3 GHz, and (b) at 8.4 GHz. top : integrated 
fliix density, middle : major axis, and bottom : position angle. The major axis increased with time from 
0.28 to 0.38 mas in Period 1. After reaching its largest value of 0.42 mas, it decreased to 0.32 mas at the last 
epoch in Period 2. The flux density also increased about 50 %, and the PA changed from 80 to 50 degrees 
in Period 1, and a possible increase in the major axis and flux density can also be seen at 2.3 GHz. 
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Fig. 12. (a) Evolution of /3app as a function of distance from the core, (b) Constraints on the viewing 
angle and speed of a jet from /3app in Period 1, according to the relativistic jet model. The horizontal dotted 
line indicates a jet speed of 0.77 c, and the vertical dotted line indicates a viewing angle of 5.3 degrees. 
Assuming that the viewing angle is always 5.3 degrees, we found that the jet accelerates from 0.77 c to 0.92 
c. Assuming that the jet speed is always 0.77 c, we found that the viewing angle changes from 5.3 to 20.8 
degrees. 
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